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• Abstract. 

We study the transition probability and coherence of a two-site system, interacting 
Q ' with an oscillator. Both properties depend on the initial preparation. The oscillator is 

prepared in a thermal state and, even though it cannot be considered as an extended 
bath, it produces decoherence because of the large number of states involved in the 
C dynamics. In the case in which the oscillator is initially displaced a coherent dynamics 

i-^J . of change entangled with oscillator modes takes place. Coherency is however degraded 

Sh I as far as the oscillator mass increases producing a increasingly large recoherence 

time. Calculations are carried on by exact diagonalization and compared with two 
semiclassical approximations. The role of the quantum effects are highlighted in the 
' long-time dynamics, where semiclassical approaches give rise to a dissipative behaviour. 

Moreover, we find that the oscillator dynamics has to be taken into account, even in a 
^sO \ semiclassical approximation, in order to reproduce a thermally activated enhancement 

■ of the transition probability. 
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& ' 1. Introduction 

The quantum dynamics of a charge, moving between two potential minima, is strongly 
influenced by the variations of the surrounding geometrical configuration. The potential 
the charge is put in, is usually produced by heavy degrees of freedom (a lattice, a 
molecular structure or an environment) evolving as well. In many cases, only one normal 
mode of the heavy system is expected to be coupled with the tunnelling charge. This 
occurs for very simple molecular or mesoscopic structure, or when the time-scales of 
heavy and light systems' dynamics are so different to allow the coupling with a single 
collective mode. 

As an example, we can consider the transport of a charge between localized sites, 
in a crystal, affected by coupling with optical phonon modes, hereafter we shall refer 
to this picture also for the choice of the notation. Moreover, there are a lot of other 
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cases formally similar to it. Another example is the electron transfer reaction, in which 
an electron is excited toward an excited molecular state which is strongly coupled to 
ionic motion. In this case, coupling can reduce the tunnelling frequency of the electron 
between the two molecular states. The result is a freezing of the electron into a definite 
excited state, in which electronic and associated ionic state are entangled pQ. As a third 
example, we may consider a single-molecular conductor made by carbon nanotube [2]. 
Here, the negative differential conductance behavior is associated with phonon-mediated 
electron tunnelling processes [3]. 

In all these cases, the behaviour of the systems qualitatively changes, as the 
temperature of the oscillator increases, giving rise to a charge transfer process which 
continuously changes from coherent quantum tunnelling to incoherent classical hopping. 
Coherence properties of tunnelling systems are now accessible to a wide class of 
experiments. By broadband absorption spectroscopy, which is able to access time- 
resolved kinetics, it is possible to detect coherent oscillations in excited-state electron 
transfer of fluorinated benzenesjl]. In a single quantum dot, Rabi oscillations have been 
detected using quantum wave function interferometry [5] . Here, the the electromagnetic 
field strongly couples with excitonic levels of the dot. 

Because of the large differences among all these physical cases, the corresponding 
tunnelling system lives in very different regions of parameters space. For example, 
in molecular systems, the oscillator frequency can be much larger than the tunnelling 
amplitude of the system, leading to an antiadiabatic regime. While, in solid state physics 
dispersionless oscillators can modelize optical phonons of the systems whose frequencies 
are often much lower than hopping amplitudes of the itinerant electrons, leading to 
adiabatic regime. 

To perform an extensive study of coherence and tunnelling, as system's parameters 
span all the accessible phase diagram, we choose a simple model of a single tunnelling 
system, coupled to a single oscillator. Such a model has been widely analyzed, in 
different regimes and approximations, because of its importance both as a building 
block of cluster expansion in a lattice model[6], and as a model for chemical reactions 
and charge transfer in organics. Following a block diagonalization technique, introduced 
by Fulton and Gouterman [7], we calculate quantum propagators. These results has 
been used for calculating exactly the finite temperature spectral functions [8] and to 
characterize the band-width behaviour with temperature. Here, we study the quantum 
dynamics of the system, considering the transition probability and coherence. For this 
purpose, we study the reduced density matrix, taking into account two different initial 
system preparations: i) the electron preparation, in which oscillator is taken from an 
thermal equilibrium distribution in absence of interaction, and ii)polaron preparation, in 
which oscillator's initial state is taken from a thermal equilibrium distribution in absence 
of tunnelling. Notice that the two considered preparations are also representative of 
systems in which the charge is promoted to a given energy level with (polaron) or 
without (electron) vibronic relaxation. Reduction is obtained by tracing out the bosonic 
degrees of freedom. The particle transitions are characterized by the diagonal elements, 
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while the degree of coherence is given by the so-called purity which give a quantitative 
indication of how the system's state is a pure quantum state. The temporal behaviour of 
these quantities depends on the coupling strenght, as well as the adiabaticity parameter, 

1. e. the ratio between tunnelling and oscillator characteristic times. 

In the antiadiabatic regime, when the oscillator is faster than the tunnelling system, 
polaron preparation guarantees a coherent behaviour up to temperatures of the order 
of the phonon frequency. On the contrary, the electron preparation gives rise to a fast 
decoherence because of the entanglement with the oscillator mode, which produces a 
dissipative effect even at zero temperature. As temperature increases, both preparations 
gives the same incoherent dynamics. 

In the adiabatic regime, the phonon spectrum tends to a continuum and polaron 
recoherence times becomes longer and longer. We observe a decoherence in both electron 
and polaron preparation. In this regime, we compare the exact result with a static (SA) 
and a Quantum-Classical (QC) approximation yet introduced in [9]. We find that QC 
approximation is able to capture the high temperature polaron incoherent motion, while 
SA is sufficiently accurate only for short times. 

As a drawback of the finiteness of the system, the equilibrium is never really reached 
and in principle infinite recoherences appear. This purely quantum behaviour is not 
recovered neither by SA nor by QC semiclassical approximations. To observe a real 
dissipation it necessary to introduce a reservoir with a continuum spectral density. 
Nevertheless, we expect that in an intermediate time-scale, between the initial dynamics, 
driven by the fast degree of freedom, and the recoherence times, our single-oscillator 
model reproduces the many-oscillator case, providing that a mode with dominant 
interaction can be separated from the rest of the bath. 

The model is described in Sec. El In Sec. |3]we introduce the reduced density matrix 
for the polaron and the electron. In Sec. 15.11 we describe the exact mapping, by means 
of Fulton- Gouterman transformations, from the original electron-phonon problem into 
two single anharmonic oscillators. Then the QC and SA approximation are described. 
In Sec. [5] are presented the results and the comparison between these three different 
techniques is discussed. Sec. [6] is devoted to the conclusion. 

2. The Model 

The model we shall consider is described by the following Hamiltonian 



describing a spin-1/2 interacting with a harmonic oscillator of frequency ujq- The model 
can be associated to a large number of physical systems [IQIQT] but, for sake of clearness, 
we shall refer to an electron, in the tight binding approximation, moving in a two-site 
lattice and interacting with it by the local distortion of the lattice site [5j. In particular, 
it can be shown that it is equivalent to the Holstein two-site model [T2l [131 [8], with 
operators a and a) referring to the relative phonon coordinate and providing that 




(1) 
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fermionic operators are mapped into a pseudo-spin notation a z = c\c\ — c 2 C2 and 
o~ x = c\c 2 + c\c\. The center of mass coordinate can easily be decoupled (for a more 
detailed discussion see (HJ[H]). Therefore throughout this paper we shall refer the word 
electron to the tunneling system and the word phonon to the oscillator. 

The strength of the electron-phonon interaction is given by the constant g — gj y/2, 
J is the electron wave function overlap or hopping and 2J is the tight binding half- 
bandwidth. 

Beside the temperature, we can choose two parameters that characterize the model 
i) the bare e-ph coupling constant A = g 2 / (u Q J) given by the ratio of the polaron energy 
(E p = —g 2 /uo) to the hopping J and ii) the adiabatic ratio 7 = too/ J. 

In terms of these parameters we can define a weak-coupling A < 1 and strong 
coupling A > 1 regimes, as well as an adiabatic 7 < 1 or anti-adiabatic 7 > 1 regimes. 

Notice that instead of choosing A as coupling constant we may choose another 
combination which is more appropriate in the so called atomic (J = 0) limit i.e. 
a = ^J\/ (27) (see Appendix). 

3. Reduced Density matrix 

The study of the charge dynamics is not trivial because, in general, it is entangled 
with the harmonic oscillator. The time dependent correlation functions of the two-site 
Holstein model has been investigated in the past [15] and also a short time transfer 
dynamics has been introduced in [16J. 

In this paper, we introduce a density matrix approach for the charge dynamics 
over a very large time range. Hereafter, we shall assume that charge and oscillator are 
initially separated, being the former localized on the first site and the latter in a mixed 
thermal state. The corresponding density matrix is 

p(o) = E— ^H^X^I®! 1 )^ (2) 

where we used the notation |1) = c{|0) and (3 is the inverse temperature. The state 
\<p n ) depends on the choice of the initial preparation [17], in this paper we study two 
different situations obtained from two different limiting regimes: 

(i) electronic preparation (el): the electron is initially free (g = 0) and the oscillator is 
at its thermal equilibrium 

P {d) (0) = J2^—\n)(n\®\l)(l\, (3) 

n 6 

(ii) polaronic preparation (pol): electron is initially localized (J = 0) on a given site 
(say 1), while the oscillator is displaced accordingly (see Appendix) 

P (po °(o) = E^^K)(^l®|i)(i|, (4) 
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The dynamics is obtained by switching on g, in the first case, and J, in the second one, 
and letting evolve the density matrix with the Hamiltonian ([1]) p(t) = e~ tHt p(0)e lHt . 
The temperature enters only in the initial state trough the incoherent distribution of 
the initial oscillator states in both preparation. 

Tracing over the oscillator degree of freedom, we obtain the electron reduced density 
matrix 

pW(t) = Tr ph {p(t)} : (5) 
which, in terms of the oscillator's number states, is 

-f3ui n 

p W(t) = -^(m\e- lHt \n, l)(n, l\e lHt \m), (6) 

n,m 

To characterize the motion of the polaron we cannot reduce the density matrix by 
tracing out the phonon degrees of freedom, this is because the polaron itself contains 
phonons. In order to understand better the polaron dynamics, let us first apply a Lang- 
Firsov transformation D (see Appendix), the new fermionic particle corresponds to a 
polaron, so the density matrix with the initial localized polaron can be written as 

p^\t) = Tr ph {&p(t)D}, (7) 

and reads, in terms of the oscillator's number states, as 

-flu) on 

~~z 



p ( P oi)^ = ^ ^— - — <7n|e-*- ff *In, l)(n, l\e iRt \m). (8) 



3.1. Quantities of interest 

In this paper, we will study two measures: one for transition probability and the other 
for the degree of coherence. The diagonal elements of the reduced density matrix, in 
the site basis, represent the population of each site. The transition probability from site 
1 to site 2 is given by 

w lj2 (t) = (2|p(t)|2). (9) 

where p is the reduced density matrix in any of the previously introduced preparations. 

The off-diagonal elements of the reduced density matrix represent the quantum 
interference between localized amplitudes. However their knowledge is not sufficient to 
determine whether the state is pure or not. Suppose that initial state is pure, if diagonal 
elements do not evolve in time, the suppression of the off-diagonal elements implies the 
evolution into a mixed state. In this particular case, the knowledge of off-diagonal 
elements also determines the purity of the system. In the more general case in which 
all the elements of p evolve, the choice of the off-diagonal elements obviously depends 
on the basis. A basis independent measure for purity (called purity itself) is 

P(t) = Ttp 2 (t). (10) 

where again p is the reduced density matrix. It is easy to see that 1/2 < P < 1 with 
P = 1 if and only if the state is pure and P = 1/2 when the state is maximally mixed. 
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The behaviour of our finite system results as a superposition of oscillations with 
many different characteristic frequencies. To disentangle the relevant timescales at a 
given time t it is found useful to consider the time averaged transition probability and 
coherence, defined as 

Q(t) = - fdt>Q{t>), (11) 
t Jo 

where Q can be either u>i,2 or P. 



4. Methods 



In this section we present the methods which we use to get the reduced density matrices 
for both initial preparations. 



4-1. Exact diagonalization 

As shown by Fulton and Gouterman [7j, a two- level system coupled to an oscillator 
in such a manner that the total Hamiltonian displays a reflection symmetry, may be 
subjected to a unitary transformation which diagonalizes the system with respect to 
the two-level subsystem [3, [18j EEl EQ] . This method can be generalized to the N-site 
situation, if the symmetry of the system is governed by an Abelian group [T9] . 

In particular, an analytic method for calculating the Green functions of the two- site 
Holstein model is given in jH [21]. Here, the Hamiltonian is diagonalized in the fermion 
subspace by applying a Fulton Gouterman (FG) transformation. So the initial problem 
is mapped into an effective anharmonic oscillator model. It is possible to introduce 
different FG transformations for the electron and the polaron. The new problem results 
to be very simplified and very suitable to be numerically implemented. Analytical 
continued-fraction results exist for the electron case[8| |2T]. 

In this section, we briefly recall the FG transformations method. The density 
matrix elements are given explicitly in terms of effective Hamiltonians and calculated 
by means of exact diagonalization 

The FG transformation we use for the electronic case is 

it a \ 

(12) 





the new Hamiltonian H = VHV^ 1 becomes diagonal in the electron subspace 

*.)• (13) 

the diagonal elements, corresponding to the bonding and antibonding sectors of the 
electron subspace, being two purely phononic Hamiltonians 

H± = uootaT J(-lf a -g((J + a). (14) 

The operator (— l) ata is the reflection operator in the vibrational subspace and it satisfies 
the condition (— l) ata a(— l) ata = —a. A wide study of the eigenvalue problem was carried 
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out in [22], both numerically and analytically, by a variational method, extending the 
former results given in [T3]. In [22] H± is approximately diagonalized by applying a 
displacement, the dynamics is reconstructed by the calculated eigenvectors and energies. 

The evaluation of the polaron Green function can be done on the same footings, 
but the expression involves also the non diagonal elements of the resolvent operators, 
causing an exponential increasing of the numerical calculations. 

To avoid this problem, we first perform the LF transformation and then apply, on 
the resulting Hamiltonian (T4l) . a different FG transformation 

1 / -i / -i \ n J n \ 



Vl V2\ (-l) ata 1 
The new Hamiltonian Hlf = ViHV{~ is 

( H, \ 



H 



LF 



H 



(15) 



(16) 



where 

H± = U a^a + J(_l)<^ e T2a(at-a) + ^ 

is real and symmetric but not tridiagonal in the basis of the harmonic oscillator, the 
matrix elements of H± are given in [8]. 

In order to write down the density matrix elements, let us introduce the following 
notation: 

(m\e- iH±t \n) 



R {±) (t) - 

M™{ n (t) 
M 2 m f(t) 



(m\e- iH±t \n) 



1 



(m, l\e- lHt \n, 1) = - [rM (t) + R^ n (t) 



(m,2\e- iHt \n,l) 



(-1)' 
2 

1 



-iHt 



(Vi, l\e~* m \^ n , 1} = - + (-ir +n Rtf n (t) 



(18) 
(19) 

(20) 
(21) 

(22) 



hi i) = 2 [(- l ) nR U( t ) - c-irsftiw] ■ ( 23 ) 



The reduced electron density matrix elements are 



pi?(t)=E 



-f3u)Qn 



\N?i n {t)f 



f&(t) = E 



(24) 



the calculation for the polaron case gives 



ff\t) = e 



\M?{ n {t)\' 



P&\t) = E 



-M™i n {t)Ml^ n {t). 



(25) 
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A qualitative insight into the relevant timescales involved in the evolutions of 
and gained by looking at the behaviour of the spectral functions of the 

model (CD) studied in our previous work [8]. In terms of the Fourier transform of function 
N™{ n (t), the electron spectral function A(ou) can be defined as 

A( U ) = _±j m £ l_JV£f("). (26) 

Analogous equation holds for the polaron spectral function relating it to the function 
M hl {u). 

An example of A(u) is reported in figure [U We notice that three energy scales 
(depicted schematically in figured]) can be associated A(uS) [8J. One is the separation 
of the low lying energy level AE, the other is the phonon energy uq and finally there 
is the tunnelling J. They are depicted schematically in figure [D These energy scales 
define three different timescales: 

i) tj = 2tcJ-\ 

ii) t U0 = 2nu 1 

iii) tq = 2-kAE" 1 . 

As it is reasonable from the relation between spectral functions and reduced density 
matrix ( l26j) . these characteristic timescales are recovered in the reduced density matrix 
evolution. 



10 +1 



10- 1 
I 10 " 2 

< 10 3 
10 4 

-4-3-2-1 1 2 3 

(!) 

Figure 1. Electron spectral function for 7 = 0.1 A = 2and T = see ref. [8]. 



4-2. The static approximation 

The case, in which a light quantum particle interacts with much more massive particles, 
is very common in solid state and molecular physics. We discuss the adiabatic regime, 
meaning that, in a characteristic time for the the light particle dynamics, the heavy 
degrees of freedom can be considered approximately quiet. Here, we describe the SA 
approach in its basic formulation for the dynamics. 
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The Hamiltonian ([T]) can be written in the coordinate-momentum representation 
p 2 mul 2 g 

H = h — — - x -=xo z - Ja x - uj , (27) 

Zm 2 y 2 

with g = g^2mujQ. In the adiabatic limit (7 1), the phonon is much slower than the 
electron (heavy phonon and large electron tunnelling amplitude) and one can neglect the 
phonon kinetic term in ( 1271) . This is the well known Born-Oppenheimer approximation. 
In practice, it consists in studying the electronic problem with x as a classical parameter. 

Within this approximation, we put loq = ^Jk/m — > (m is the ion mass) and the 
Hamiltonian becomes 

Had = ^ ~ 7/| xcr z - J&x- (28) 

The eigenvalues can be expressed trough the classical displacement x 

V ± (x) = ^x 2 ±Q(x), (29) 



with Q(x) = \j\x 2 + J 2 . The lowest branch (— ) of ff29|) defines an adiabatic potential 
which has a minimum at x = as far as A < 1 while for A > 1, it becomes double well 
potential with minima at ix m , x m — \J ^7 — ^ case the electron is mostly 
localized on a given site. The quantum fluctuations are able to restore the symmetry in 
analogy to what happens for an infinite lattice [23J. It is worth noticing that, in this limit, 
Hamiltonian ([1]) is equivalent to the adiabatic version of the spin-boson Hamiltonian 

[SUES]. 

The temporal evolution is given by 

e -iH ad t = e -i*£*[co S n(a;)t + i( ^ a z + -^—a x ) smn(x)t], (30) 

V2i2(a;) \l{x) 

so the density matrix dynamics can be explicitly calculated 

The electronic initial preparation, corresponds to the density matrix 

p(0) = |1)(1U/^ / dxe-^\x)(x\, (31) 
V 2n J 

tracing out the phonon we obtain the electron reduced density matrix with elements 



(5J\ r , _ pjxo sin 2 (Jt J u 2 X 2 + 1 

/ due. 2 1 



P2S = I due 



1 + \ 2 u 2 



p[e? = ,pjx f due _^HWu 2 x 2 + i) (32) 



2tt J 2x/A 



2 u 2 



where the scaled lenght u = xk\/2/g was introduced. 

In the same way, we can introduce the polaronic preparation 



p(0) = |l)(l| e -ift-4/^ / dxe~^ x2 ~^ x) \x)(x\, (33) 




It is worth noting that, in the adiabatic limit, we cannot define the polaronic dynamics, 
as introduced in (J7J), because the operator D is not defined for luq = 0. In this limit, 
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we study electronic dynamics with an initial polaronic preparation. The corresponding 
reduced density matrix is 



P2,2 



JpoO 
r 1,2 



I/3JX 
~2^~ 



du e 



-^(«-i) 2 ' 



due-^^ 2 



sm 2 (Jt^/(uX) 2 + 1) 
{uX) 2 + 1 
uX sin 2 (Jt J '{Xu) 2 



((Xuy + 1] 



— %- 



sin(2JtJ(Xu) 2 + 1) 



2J(Xu) 2 + 1 



(34) 



It is possible to show that P22 is actually the adiabatic limit of the diagonal 
element of the reduced polaronic density matrix, while this is not true for the off- 
diagonal elements. 

We want to stress that, in the SA approach, the phonon is completely static because 
its momentum p has been neglected. Here, only the initial phonon distribution plays a 
role, but during electron hopping, oscillator is taken to be fixed. 



4-3. A quantum- classical dynamics approximation 

To account for dynamics of the slow variable, a mixed quantum-classical dynamics 
can be introduced. In the past, several schemes for quantum-classical dynamics 
has been proposed, for example starting from the Born-Oppenheimer (SA) adiabatic 
approximation for the ground state at each step and using a density functional 
Hamiltonian [261 I2Z]- Another approach, good for a short time dynamics, consists 
in a mapping from the Heisenberg equations to a classical evolution by an average over 
the initial condition [28l [29] . Some schemes are based on the evolution of the density 
matrix coupled to a classical bath [301 IS]- A systematic expansion over the mass ratio 
has also been done, starting from partial Wigner transform of the Liouville operator, in 
[32l [33l [34] . The QC approximation we use is essentially that of refs. [30~1 l3T]. 

Let us consider Hamiltonian ( 1271) . where x and p are assumed to be classical 
variables which can be represented as the components of a vector u. Then a QC state 
vector can be introduced as 

/ x \ 



u <g) a 



x 

P 
o- x 



(35) 



The classical variables evolve with the Ehrenfest equations 
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while the quantum variables evolves in the Heisenberg picture 
& x = -V2gxa y 



& y = V2gxa x - 2Ja z . 



(37) 



2Ja„ 



To give a unified description of the overall evolution, we define a Liouvillian operator 



£ = C x + C p + L a with 



( -V2gx \ 

-i V2gx -2J (38) 

v 2J j 

and £ x = x-^ £ p = p-j^. So, the time evolution is given by 

v(t) = e l£ 'v(0). (39) 

The numerical integration can be implemented using the symmetries Trotter 
breakup formula [351 EB] 



N 



v(0) 



(40) 



v(£) ~ (^ e iC ' T ^e iCp 2e iCxt e iCp 2e i 

with e = t/AT. All the density matrix elements, can be expressed in terms of elements 
of v(t). 



5. Results 



5.1. Antiadiabatic regime 

In figure[2]is shown the time behaviour of the purity P (Eq. (|10p ) as well as the transition 
probability ( (jHj)) obtained in the antiadiabatic regime when the phonon frequency (o;o) 
is much larger than electron hopping J for both (el) ([3]) and (pol) (j3j) initial preparations. 
Timescales defined in section are shown as vertical lines, the time scale is logarithmic 
to better show the much different time domains. We consider two parameters sets at 
several temperatures. One characteristic of strong coupling (left panels) and the other 
of weak coupling (right panel) . The same sets of parameters and temperatures is used 
in figure [3] where with show the time-averaged P and w. Let us first discuss the strong 
coupling regime. 

It is known that, in the antiadiabatic regime, the polaron is a well defined quasi- 
particle at strong coupling [37], in the sense that, in the polaronic spectral function, 
almost all the spectral weight is contained in the polaronic peak. This has also been 
shown for a two-site model [8j [HI [15], [381 ESS SO]- On the contrary, in the electron 
spectral function, the total spectral weight is distributed between a large number of 
frequencies [S]. 

From the point of view of transition probability and purity, the strong dependence 
on initial preparation can be seen comparing the low temperature evolution for both 
the actual (figure [2] left panel) and the time-averaged (figure [3] left panel) quantities. 
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Let us consider the electron preparation (Figs. EIE! bottom left panels). Here the 
coupling with the oscillator mode is strong and the system starting from a disentangled 
state (equation (jHJ)) evolves in a state in which electron is entangled with the oscillator. 
This is shown in the purity evolution where we see that electron loses coherence very 
rapidly, on a time scale tj, and becomes a mixed ensemble, even at zero temperature. |jj 

On the same timescale, transition probability approaches 1/2 in average (figure [3] 
bottom left panel). The initial decoherence is almost independent on the temperature, as 
can be seen from averaged quantities, while long time recoherence peaks are suppressed 
as T increases. Such a suppression results from destructive interference between 
time evolution of the different terms appearing in (J6J) when excited oscillator states 
are initially populated. Referring to the spectral analysis [8] and to figure [Q, this 
phenomenon must be ascribed to the superposition of a large number of high frequency 
excitations. 

Even if the transfer does not have a regular shape, one can see some high frequency 
oscillations of period t Uo . These frequencies correspond to the energy separation between 
two adjacent electronic bands [8|. 

On the contrary, the polaron preparation (see figure El top left panel) evolves in 
a state which is completely coherent at zero temperature. The frequency associated to 
polaron transfer is equal to the renormalized tunnelling J*, as predicted by the HLFA 
(see [6]). So, the state is pure and delocalized. The polaron state remains coherent 
even for temperatures comparable with uj , but higher frequency modulation appears 
making the state oscillate from a pure to a mixed one. Nevertheless, it is possible to see 
an overall modulation of the transition probability with the same period tq even at the 
largest temperature. This is in contrast with the HLFA at T ^ (CSJ) which predicts that 
the polaron band decreases with temperature and consequently tq increases. However 
the purity decreases as temperature increases, as shown in figure El This is an effect 
of the broadening of the polaron band that is observed as temperature increases [8]. 
Indeed, a distribution of spectral weight among several poles around the polaron band 
occurs as an effect of increasing vibronic excitations (ref. [8] figure 3 upper panel). This 
leads to a decoherence effect due to destructive interference between these oscillating 
contributions to purity (see Eqs. (118f25p ). For high temperature (T ^> uq), the state 
becomes completely mixed and the evolution of the polaron is analogous to that of the 
electron. This is evident from the highest temperatures curves shown in figure [31 left 
panels upper and bottom left. We conclude that the main source of decoherence is 
temperature for polaron, while the electron decoheres even at zero temperature due to 
the coupling with the vibronic mode. 

This is also found in the weak coupling regime (electron preparation Figs. I2ll3l 
bottom right panels). Here electron coherence approaches a value which is larger in 

| Notice that, in this case, the analysis of purity and transition probability alone in principle do not 
allow to determine which states the mixture is composed of (pointer states). In particular whether the 
states are localized or not. However, a straightforward analysis of non diagonal elements of reduced 
density matrix shows that the states are indeed localized. 
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Figure 2. Polaron (top) and electron (bottom) populations and purity. Left panels 
antiadiabatic and strong coupling regime: 7 = 10 A = 40. Right panels antiadiabatic 
and weak coupling regime 7 = 10 A = 10. Curves are for T/ujq = 0.0 (black), 
T/ujq = 0.5 (blue), T/u = 2.0 (green), T/oj = 10.0 (red). Vertical lines marks 
from left to right the timescales t Uo ,tj,tq. 

average than that at strong coupling (figure [3]) and decreases as temperature increases. 
We see (Figs. [2|3l right panels) that the polaron purity differs qualitatively from the 
electron only near zero temperature, while averaged polaron and electron transition 
probabilities is essentially the same for all showed temperatures. Indeed, a sufficiently 
weak interaction is not able, at zero temperature, to excite many vibrational states, 
so the electron decoherence is essentially given by the small perturbation of the lowest 
oscillator's states. On the other hand, the polaron is not formed (we are below the 
polaron crossover) and the charge does not acquire much coherence by moving with 
the oscillation cloud. Spectral analysis shows (ref. [8] figure 3 lower panel) that in 
weak coupling HLFA is qualitatively recovered, we have a polaron band narrowing as 
temperature is raised up in contrast with the strong coupling behaviour. 

It is worth stressing that, in both weak and strong coupling regimes, at high 
temperature, the increasing number of oscillator states involved in the initial state 
produce decoherence on timescale tj. Decoherence can be partial but nonetheless 
no environment is needed to explain the decoherence process. The only source of 
decoherence are the states populated by the initial thermal distribution 
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Figure 3. Polaron (top) and electron (bottom) time-averaged populations and purity. 
Parameters and labels are the same of figure [5] 

5.2. Adiabatic regime 

Results form the Exact Diagonalization method is reported in figure HI the averaged 
quantities are shown in [5] in the same way we did in the antiadiabatic case. Notice that 
now the shortest timescale is tj. 

Let us first discuss the strong coupling case. We see that, in contrast with 
antiadiabatic regime, there is a marked dependency on temperature of both electron and 
polaron properties. More specifically, polaron preparation no longer evolve coherently 
at low temperature. 

In the first time scale, Tj, the particle is localized (its transition probability is 
extremely low) but the state keeps on being quite pure. The polaron is trapped 
inside the initial site and both transition probability and coherence evolve initially with 
characteristic time tj independently on temperature. At intermediate timescale t uo 
temperature induces derealization while the coherence decreases. In this time regime, 
the polaron transition probability is related to the quasi classical motion of the oscillator 
and depends strongly on temperature. 

This can be seen in figure El where we plot the temperature dependence of the level 
reached by the averaged transition probability on the time scale t Wo (inset). Since there 
is no clear plateau in the averaged transition probability for times greater than t Uo , the 
choice of the transition probability level it is rather arbitrary. We choose the value of 
W\2 at r^ . We see that this transition probability level pass from a low temperature 




Figure 4. Polaron (top) and electron (bottom) populations and purity. Left panels 
Adiabatic and strong coupling regime: 7 = 0.1 A = 2. Right panels adiabatic and 
weak coupling regime 7 = 0.1 A = 0.5. Curves are for T/uj q — 0.0 (black), T/ui = 0.5 
(blue), T/u)q = 2.0 (green), T/ujq = 10.0 (red) .Vertical lines marks from left to right 
the timescales tj,t Uq ,tq. 



behaviour, which is temperature independent, to a temperature dependent behaviour 
trough a wide crossover. 

At very low temperature, after the characteristic time t Wo , coherence and transition 
probability reach a quasi stationary value that is essentially dominated by fast tunnelling 
of the charge between the two sites, with a given phonon displacement. Once 
temperature increases, classical activation processes of the phonon coordinate becomes 
effective, producing an increase in transition probability as well as a decrease of the 
purity As we shall see in the next section, this thermal activated behaviour is to be 
ascribed to the phonon classical hopping between two adiabatic minima and disappears 
in the SA where such hopping events are absent. 

As far as electron is concerned, we can see that the site occupancy begins to 
oscillate coherently, with period tj, with a damping increasing with temperature. In 
the same timescale, the time averaged value shows a saturation at low temperature. At 
a temperature independent intermediate timescale r^ , time averaged coherence reaches 
a very slowly decreasing level which decreases with increasing temperature. Quantum 
oscillations still exist further in time, but the remaining coherence slows down at high 
temperature, so that the long timescale (tq) seems to be not relevant in this case. The 
electron's tendency to coherently hop is suppressed by decoherence, induced by excited 
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Figure 5. Polaron (top) and electron (bottom) time-averaged populations and purity. 
Parameters and labels are the same of figure 2] 




Figure 6. Time averaged transition probability for A = 2.0,7 — 0.1, temperatures are 
from bottom to upper curves T/ujq = 0.1,0.5, 1.0,2.0, 10.0,20.0. Vertical lines marks 
from left to right timescales tj,t Wo and tq respectively. Inset: levels reached at time 
t~ui (arrows in the main panel) as a function of the inverse temperature. 
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phonons, whose number increases with temperature. In the adiabatic limit, this effect is 
more evident than in the antiadiabatic case, this is because the energy spacing between 
the oscillator's levels becomes very small and the spectrum tends to a continuum. 

In the weak coupling case (figure H] and [5] right panels) the transition probability is 
very similar for both electron and polaron as in the antiadiabatic case. Conversely, at 
low temperature, the electron is much coherent than the polaron. In this regime, being 
the adiabatic potential single-well, the displaced phonon base is not the best choice. So, 
many displaced oscillator states are coupled with polaron, which decoheres rapidly. 

5.3. Comparison with quanto- classical approaches 




Figure 7. Polaron (top) and electron (bottom) populations and purity in the adiabatic 
strong coupling regime 7 = 0.1 A = 2. Left panels low temperature T/u> — 0.1. Right 
panels high temperatures T/luq = 10.0. Curves refers to ED (black), QC (blue), SA 
(green) approximations. Vertical lines marks from left to right the timescales tj ^t^^tq. 

In this section, we show a comparison between the results obtained in the three 
different ways described before: exact diagonalization (ED) by means of mapping 
introduced in Sec. I5.1[ the quantum-classical (QC) dynamics approach described in 
Sect. 14.31 and the static (SA) approximation (Sect. I4.2p . We shall limit ourselves to an 
adiabatic case (7 = 0.1) with electron-phonon interaction strong enough to allow the 
polaron formation (A = 2). 

In figure d is reported the exact dynamics given by the three different techniques, 
while in figure M is shown the time average. Remember that, as far as the polaron is 
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Figure 8. Polaron (top) and electron (bottom) time-averaged populations and purity. 
Parameters and labels are the same of figure [7] 
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Figure 9. Time average transition probability level at the timescale t Uo in the 
adiabatic regime (A = 2.0,7 = 0-1) as a function of inverse temperature for different 
preparations and different approximations. 

concerned, in both SA and QC approximation, the purity represents that of the electron 
with a initially displaced phonon distribution. 

At low temperature (left panels), and within the tj timescale, the classical phonon 
is almost freezed, and so both SA and QC approximation are equivalent. Nevertheless, 
the ED behaviour, of both polaron and electron preparation, is quite different because 
of quantum fluctuations. In particular, for short timescales, one can see that the ED 
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dynamics, at low temperature, is damped faster than the other two approximations. The 
difference becomes more evident for higher timescales. In this regime, the temperature 
is not so effective in dissipation processes, while the strong coupling and the quantum 
uncertainty produce a sort of purely quantum damping. 

A damping of this sort was also found in ref. p2], where authors consider a spin- 
boson Hamiltonian for a tunneling system coupled to a multi-mode bath with a ohmic 
spectral density. It is worth stressing that, in the present case, such a damping is not 
due to the bath but rather a consequence of the entanglement between the tunnelling 
system and the quantum oscillator. Nevertheless our results are comparable with that 
of ref. [I?]. The reason is that, at strong coupling, the charge is coupled to the bath 
through a single collective mode, taking into account all the bath's oscillators. 

As expected, the exact polaron purity completely differs from the purity obtained 
by SA and QC. The reason is that, at low temperature and strong coupling, the polaron 
is well defined and moves as a quite coherent particle, as can be seen from the long 
timescale oscillations. On the contrary, if we consider only an electron with a displaced 
oscillator, the particle remains localized because the trapping mechanism, but it cannot 
coherently tunnel between the two sites. 

A semiclassical behaviour is approached for T greater than uu , when SA and QC 
reproduce the ED transition probability within the tj timescale, as results evident in the 
right panels of Figs. (JTHHJ) - The QC approach remains a good description also for higher 
timescales. Notice that the same occurs in the transport of extended system where 
classical incoherent transport is achieved when T is greater than 0.2uo [S]. It is worth 
noticing that for very high temperature (T ~ J,T/ujq = 10 in figs. [8|9|) the polaron 
is not formed, its dynamics approaches that of the electron in an initially delocalized 
phonon distribution. As a result, the QC purity approaches the ED's. 

At high temperature, the oscillator dynamics plays a relevant role, QC is a much 
better approximation of ED than SA. This fact can be understood by realizing that the 
main temperature effect is the damping of the coherent tunnelling oscillations. Once 
these oscillation are sufficiently suppressed, the phonon driven dynamics prevails. In 
the SA framework, the initial thermal distribution of the phonon coordinate makes the 
electron thermalizes irreversibly in a time that is the shorter the greater the temperature. 
Before this adiabatic thermalization, i.e, in a tunnelling period, the SA is still a good 
approximation. 

Afterwards, the hopping of the oscillator coordinate into the other minimum of the 
adiabatic potential (equation fl29|) ) takes place. The charge degree of freedom follows 
while wi t 2 saturates on average. In extended systems this regime corresponds activated 
mobility regime pi2l Ff3l WQ . Since SA completely neglects the oscillator's dynamics, 
it does not predict correctly u>i j2 , as can be seen in figure [9j The QC approximation, 
instead, gives a correct qualitative prediction. 
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6. Conclusions 

In this paper, we have studied a simplified model to treat the dynamics of a tunnelling 
charge interacting with a vibrational degree of freedom. We introduced a reduced density 
matrix approach to characterize the charge dynamics. Temperature is introduced 
by taking an initial equilibrium distribution of the oscillator. Both the transition 
probability and the purity are studied, in order to connect the charge transfer with 
its coherence. 

Due to the simplicity of our model, we were able to span all the parameter's space 
even at high temperature and strong coupling and to study the role of initial preparation. 
Moreover, we can explore a temporal range which is very large, compared with the 
typical time scales that can be obtained in models where the charge is coupled with a 
many degree of freedom oscillators' bath [TTJ H3, H6] . 

As in any finite system, in our model, transition probability and purity can be 
expressed as a superposition of many non commensurate by oscillations. We therefore 
expect an oscillatory behaviour in our quantities of interests. However the initial thermal 
distribution of the oscillator states induces decoherence on intermediate timescale, due 
to the strong interaction with the oscillator. This phenomena occurs depending on 
initial preparation of the system. 

We find that, in the antiadiabatic and strong coupling regime, the polaron exhibits 
a coherent tunnelling dynamics over time scales of the inverse polaron renormalized 
band. The coherent behaviour is lost out of the polaronic phase, i.e. increasing the 
temperature or decreasing the coupling. Electron evolves though partially incoherent 
dynamics. In the adiabatic strong coupling regime, temperature enhances the incoherent 
polaron charge transfer. The opposite occurs in the electron preparation. 

In the adiabatic regime, two common approximations has been compared with 
exact results, the aim is to highlight the limits of validity of these approximations 
and to provide a simple testing tool, the two-site model, for generalizations to other 
extended models. As expected, a dynamical semi-classical approximation gives good 
estimates for both coherence and tunnelling amplitude, at high temperature T 3> u> . 
Quite unexpectedly, it allows for a good approximation of the transition probability 
at low temperature, as far as time averaged quantities are concerned. However, such 
a quasiclassical approximation fails approaching the anti-adiabatic regime where non 
adiabatic transition are expected to contribute significantly to charge dynamics. This 
simplified model could serve to test approximate schemes to deal with this regime |45j. 

To conclude, we have shown that a non dissipative evolution of a tunneling system, 
strongly coupled to a single oscillator, can give rise to decoherence phenomena when 
the initial distribution of the oscillator is thermal and when the oscillator distribution 
is not initially equilibrated in the presence of the charge, that is sufficiently far from 
the thermal equilibrium distribution in the presence of the charge. These decoherence 
phenomena are independent on the presence of a dissipative bath. Thus, in a non 
equilibrium experiment in which a charge is introduced in a molecular system and 
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interacts strongly with a particular mode of the molecular system, decoherence effects 
can be triggered alone by this coupling and by the initial non equilibrium distribution 
of the molecule. 

Acknowledgments 

Authors wish to thank Sara Bonella, Carlo Pierleoni and Sergio Caprara for fruitful 
discussions and suggestions. 

Appendix 
Atomic limit 

In the atomic limit, the Hamiltonian is diagonalized by the so-called Lang-Firsov (LF) 
transformation 

D = e aa ^~ a \ (.1) 

This transformation shifts the phonon operators by a quantity a, while the electron 
operators are transformed in new fermionic operators, with energy E p , associated to 
a quasi-particle called polaron[17l [12]. This particle can be tough as a charge moving 
together with a dressing cloud of oscillator quanta, a 2 represents the mean number of 
phonons in the polaron cloud. 
The atomic Hamiltonian 

H = u>ocra — go~ z {a) + a), (.2) 

after the LF transformation H = D^HqD becomes: 

H = cu a?a + E p /2, (.3) 

the eigenvalues E n = Uqti + E p /2 correspond to the two-fold degenerate eigenvectors 
iV'n'i) = D\n, j) = Cj\ n )> were the index n = 0, . . . , oo refers to the photon number, 
j = 1,2 to the electron site and cj is the polaron creation operator cj = DcjD^ = 
c'j exp{(— iya(al — a)}. 

In the case of finite J, the hopping term is not diagonalized by ([T]) and the new 
Hamiltonian H = D^HD becomes 

H = u)Q(va — J(o~ x cosh(2a(a^ — a)) + 

+ 'uT y sinh(2a(a f - a))) + E p /2. (.4) 

Depending on the choice of the parameters, the problem will be better described by 
a electron or polaron excitation picture. In particular, in the weak coupling limit, both 
the small polaron and the electron are good quasiparticle while, in the intermediate and 
strong coupling regimes, the polaron behaviour prevails [39J. 

The different regimes was widely studied, in literature, both for the two site problem 
and the extended case. The anti-adiabatic case was first studied in small J perturbation 
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regime [T2"l US] and in the Holstein-Lang-Firsov approximation [T21 EZ] (HLFA), were 
an effective Hamiltonian is introduced to eliminate the phonon states. In HLFA J is 
substituted by an effective hopping integral obtained by averaging the displacement 
exp [2a(cr — a)} on the thermal distribution of phonons. The resulting effective hopping 
integral is 



where n% is the Bose occupation number. At zero temperature, the well known 
exponential reduction of the bandwidth is obtained 



and, as the temperature is increased (T/J » 7), the bandwidth decreases rapidly. 
As we shown in ([8]) and in the present paper, this is a good approximation at zero 
temperature but it becomes inadequate at finite T where incoherent processes turns out 
to be important. 
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